Correlated dynamics of weakly charged silica spheres at an air-water interface* 
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Abstract - Optical microscopy and multi-particle tracking are used to investigate the spatially 
correlated motion ol weakly charged silica spheres at an air- water interlace lor different area fraction 
n occupied by the particles. When the area fraction is very small, e.g. n — 0.03, the correlation 
lunction along the line joining the centers ol particles D rr decays with inter-particle distance R as 
l/i?°- 86±0 ' 02 , and the lunction perpendicular to this line D 0e decays with R as l/R 1A5±om , which 
differs from the results ol [Phys. Rev. Lett. 97, 176001 (2006)] with low surface viscosity (where 
D rr oc l/R, Dee oc l/R 2 ). We argue that the differences arise from the Coulomb interaction between 
particles. The Coulomb interaction enhances the correlated motion ol particles. Experimental 
results show that with the increase ol n, the decay rate ol D rr and Dee with R decreases and the 
cross-correlation enhances lor the Coulomb interaction increases. The Coulomb interaction between 
colloidal particles should serve as an effective surface viscoelastical role in our system. With the 

/ \ 3/2 

scaled separation -S I j t the correlated motions lor various values ol n and different particles 

can be scaled onto a single master curve, where d is particles' diameter, r) w is the viscosity ol the 
water, and ?7es,2p is the effective surface viscosity whose measurements agree well with that ol one- 
particle surface viscosity T) e s,ip- The effective surface viscosity r] e s,2p as a lunction ol the area fraction 
n lor different silica spheres is presented. 

PACS numbers: 82.70.Dd, 68.05. Gh, 05.40.-a, 83.85.Jn 



I. INTRODUCTION 

Colloidal particles moving through a fluid create a flow that affects the movement of other particles in its vicinity. 
This interaction between colloidal particles is referred to as hydrodynamic interaction. The hydrodynamic interaction 
of bulk suspensions is long-ranged and decays with inter- particle distance R as 1 / R 0, 0| ■ In many real circumstances, 
particles are not completely free, but spatially confined in some special environments, such as microfluidic devices, 
porous media, fluid interface and cell membrane [||-[5j]. The hydrodynamics behaviors of confined colloidal suspensions 
are one of the current investigation focus. In 2000, E. R. Dufrcsne et al. [6j measured the dynamic coupling of 
particle pair above a single flat wall using optical tweezers. Their experimental results showed that while the distance 
of particles to the wall is about three times of particle radius, the hydrodynamics interaction between particle pair 
decays with R as l/R 2 . In 2011, P. L. Pushkar et al. Q made a further investigation on many-body hydrodynamic 
interaction in the presence of a single wall by holographic optical tweezers and dynamic simulation. They found 
that the behaviors of hydrodynamic interactions between particle pair and multi-particle clusters are similar, i.e., 
the hydrodynamics interaction between particles in clusters also decays with R as l/R 2 . This results differ from 
that in bulk in which hydrodynamic screen appears in high particle concentration. B. Cui et al. investigated the 
hydrodynamic interaction between particles in a narrow channel Q and confined between two plates (9[ in 2002 
and 2004, respectively. They found that a linear channel sharply screens the hydrodynamic interaction; and the 
hydrodynamic interaction of particles confined between two plane walls decays with R as l/R 2 . Perpendicular to the 
line joining the centers of particles confined between two plates the cross-correlation is negative, i.e., particles exert 
anti-drag on one another in this direction. The circulation current is the direct cause for this negative coupling. The 
above works revealed the corrected dynamic characters of colloidal particles confined by hard wall. 

There is a non-slip boundary condition near a hard wall which could cut off the fluid field in its vicinity. The fluid 
interface is a slip boundary and may not cut off the fluid field. The influences of boundary condition reflect directly 
on the correlated dynamics of particles. Because of the effects of the spatial symmetry breaking, the hydrodynamic 
interaction of particles at an air-water interface is different from that in symmetrical bulk. In 2006, V. Prasad et 
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al. [lOj investigated the hydrodynamic interaction between particles at an air- water interface and found that there 
is a transition from three-dimension (3D) to two-dimension (2D) hydrodynamics with increasing surface viscosity rj s . 
At high rj s , the correlated motion along the line joining the centers of particles D rr , and ones perpendicular to this 
line Dgg are nearly equal and have logarithmic R dependence. With the decrease of r] s , the decay rate of D rr and 
Dgg with R increases and the their values begin to deviate. At low r) s , D rr oc 1/R, Dgg oc 1/R 2 . Notwithstanding 
the differences in the behavior of D rr , Dgg at high and low r] s , all the data sets can be scaled onto a single master 
curve by rj s . Evidently, the surface viscosity ry s has a substantial influence on both D rr and Dgg. In 2011, M. H. 



Lee et al. 11| investigated the correlated motion of particles at an oil- water interface. Their experiment results also 



show a transformation from 3D to 2D hydrodynamics with the increase of rj s . The results of Prasad et al. have the 
effect of fluid field only from one side; and ones of Lee et al. have the influence of the fluid field from two side. The 
hydrodynamic interaction between particles at an oil- water interface is influenced by the viscosity of the oil r\ oi \ , water 
r\ w and interface rji nt . From the above, we know that the dynamic characters of colloidal system in the confined 
condition exhibited much more rich physical image than that in the bulk. The investigations on the dynamics of 
confined colloidal suspensions have significant application in many important field, such as pharmacy and Biological 
science etc. (l2| . 

We are interested in the influence of the hydrodynamic and Coulomb interactions on dynamic characters of colloidal 
system when the boundary condition changes from the symmetry to completely asymmetry, such as air-water interface. 
In the previous experiment investigations about the hydrodynamic interactions between particles at an air(oil)-water 
interface the interface viscosity rji n t was treated as an adjustment parameters. In this paper, we measured the cross 
correlation of colloidal interfacial diffusion under the different interfacial concentrations of weakly charged colloidal 
spheres, and investigate the influence of Coulomb interaction between the spheres on the correlated motion. 

Getting rid of the impurity from the aqueous surface, we make a clean air-water interface and investigate the 
correlated motions of silica spheres SI (d = 1.57/im), S2 (d = 0.97/im) and S3 (d — 0.73^m) at this pure interface by 
optical microscopy and multi-particle tracking. The experiment results show that the cross correlation of particles at 
pure air- water interface is stronger than that exhibited in the experiments of Prasad et al. [l(| in the case of low r] s . 
The correlated motion D rr cx i/i?0-86±o.02 turns into D ^ K 1 / i? o.68±o.02 when n wag from Q Q3 tQ q 35 The par ti c les 

used in our experiment are weakly charged silica spheres. Coulomb interactions between particles enhanced the cross 
correlation of particles. With the increase of particles' area fraction n, the mean particle separation decreases and 
Coulomb interactions increase. Thus the cross correlation of particles enhances and the decay rate of D rr and Dgg 
with R decreases as n is increased. The correlated motion of particles as a function of inter-particle distance R for 
different n is presented. As discussed in the paper, the Coulomb interaction between colloidal particles should serve 
as an effective surface viscoelastical role in our system. Though the differences in the behavior of correlated motion 

/ \3/2 

for different n and particles, with the scaled separation ( — ) and correlation function < D rr gg/r > /2D S , 

a y Ves ,2p / ' ' ' 

all the correlated motion can be scaled onto a single master curve. Where D s = ksT ' I^Kr\ e s,2-p and r\ es ,2p is termed 

as two-particle effective surface viscosity whose measurements agree well with that of one-particle surface viscosity 

Ves,ip- We present the effective surface viscosity rj eSt 2 P as a function of the area fraction n for different silica spheres. 

Particles S3 exhibit an anomalous long-range correlation which we do not understand at present. 

The paper is organized as follows. Section II is devoted to the description of the experiment setup and results. In 

section III, data analysis are presented. In section IV, further discussions are given. Finally, the work is summarized 

in section V. 



II. EXPERIMENT SETUP AND RESULTS 



The particles used in the experiment are silica sphere 1 (Sl)with the diameter d = 1.57 ± 0.06/ito, silica sphere 2 
(S2) with the diameter d — 0.97 ± 0.05fim and silica sphere 3 (S3) with the diameter d = 0.73 ± 0.04^m. The SI and 
S3 particles were purchased from Duke Scientific and were synthesized by the Stober method, followed by a sintering 
process at 900°C to drive off all organic materials and water. The particles S2 were purchased from Bangs Laboratories 
and are non porous spheres. We follow the same procedures as described in Ref. [13] for the experimental setup and 
the sample preparation and cleaning. The sample cell is made from a stainless steel disk (which is hydrophilic) with 
a central hole of diameter 13mm. The bottom of the hole is sealed with a 0.1mm thick glass cover slip, which also 
serves as an optical window. The sidewall of the hole together with the bottom glass slip forms a container, which has 
an effective height of 1.0mm. The top of the container has a sharp circular edge, which is used to pin the air- water 
interface and reduce unwanted drift. The entire cell is placed on the sample stage of an inverted microscope, so that 
the motion of the interfacial particles can be viewed from below and recorded by a digital camera. 

We take great care to clean the particle samples and the air-water interface. Surface pressure measurements reveal 
that no detectable impurity is found in the cleaned particle samples. The particle-methanol solution is then injected 
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onto a cleaned interface using a syringe pump. The silica spheres disperse well at the clean air-water interface. 
Individual particles undergo vigorous Brownian motion and remain stable at the interface for days. These particles 
remain in focus under high magnification, indicating that the silica spheres are strongly bound to the interface and 
their vertical position is determined by an energy minimum. The usual capillary effect is not applicable to micron- 
sized spheres for the gravitational energy of the particles is much smaller than the relevant energy of the interface, and 
particles can choose an equilibrium vertical position without introducing long-range deformations to the interface. 

To obtain good images of the silica particles phase contrast microscopy is used. Using homemade software with 
a spatial resolution of 60-100nm we construct particle trajectories from the consecutive images. From the particle 
trajectories, we determine the vector displacements of the tracers Ar(t,r) = r(t + r) — r(t), where t is the absolute 
time and r is the lag time. In the experiment, we measure the mean cross-correlated particle motion [2]: 

D aP (r,T) =< Ar^,T)Ar^,T)5(r-i^(i)) (1) 

where i,j are particle indicates, a and (3 represent different coordinates, and R l,J is the distance between particle i 
and j. The off-diagonal elements are assumed to be uncorrelated. We focus on the diagonal elements of this tensor 
product: D rr which indicates the correlated motion along the line joining the centers of particles, and Dgg which 
represents the correlated motion perpendicular to this line. 

The correlated motion D rr and Dgg are the function of R and r. The lag time r need approach to zero to eliminate 
its effect on R = R(r). There are two ways to calculate the r-independent quantities < D rr /r > T and < Dgg/r > r , 
which depend only on R. Firstly, if the data shows D rr , Dgg oc r within the observation period of r, the slopes of 
D rr , Dgg versus r are given by D rr /r and Dgg/r. Otherwise, we calculate the partial differential of dD rr /dr and 
dDgg/dr from which we derived the corresponding slope value at r = by the extrapolation. 

We measured the curve of D rr /r and Dgg/r vs. R for dozens of different area fraction in range of 0.01-0.35 for the 
data of SI. We make a coarse grain process by grouping the dozen of curves to range: 0.01-0.05, 0.05-0.1, 0.1-0.14, 
0.15-0.17, 0.2, 0.35 (see the figures below). We average curves in each groups respectively. The similar processing 
procedure had also been down for the data of S2 and S3. The results are obtained by averaging over 10 6 particles, 
ensuring that the statistical averaging is adequate. 

Figure 1 exhibits the measured values of < D rr /r > and < Dgg/r > for silica microspheres SI scaled by Do = 
kBT/3iTr] w d as a function of R/d with the area fraction n ranging from 0.01 to 0.05, where r\ w is the bulk viscosity of 
water in room temperature. The average result of these curves were drawn in the inset of Fig. 1 , while the corresponding 
area fraction n was estimated as 0.03. Solid and open symbols represent < D rr /r > /D and < Dgg/r > /D , 
respectively. The error bar in the inset of Fig. 1 indicates only the error range of the values of correlated motion, but 
it's slope. With the increase of inter-particle distance R, < D rr /r > and < Dgg/r > decrease following the power 
law. As shown in Fig. h< D rr /r >oc i/ij0.86±o.02 an( j < £) gg / T >oc j_^i.45±o.03 w hj c h differs from experiment 

results of Prasad et al. [10( with low surface viscosity where D rr oc 1/R, Dgg oc 1/R 2 . 

Figure 2 (a) exhibits the measured values of < D rr /r > and < Dgg/r > scaled by Do as a function of R/d 
with various estimated values of area fraction n for silica microspheres SI. Different scattering symbol represents 
the data with different area fraction of particles. As shown in Fig. 2 (a), with the increase of area fraction n, the 
magnitude of curves of < D rr /r > and < Dgg/r > increase. We fit the data in Fig. 2(a) by < D rr /r >oc l/R Xr and 

< Dgg/r >oc l/R Xe . Figure 2 (b) shows the fitted parameters A r and Xg, which characters the decay rate of cross 
correlation between particles, as a function of n. As shown in Fig. 2 (b), with the increase of n, X r and Xg decrease 
indicating that the decay rate of correlated motion vs. R decreases. Thus, the cross correlation of colloidal diffusion 
was enhanced with increasing area fraction n. 

Figure 3 (a) and (b) shows the measured values of < D rr /r > /Do and < Dgg/r > /Dq vs. R/d with various values 
of area fraction n for silica microspheres S2 and S3, respectively. The data processing for S2 and S3 is similar to that 
of SI. Every curve is derived by averaging the data with adjacent area fraction. The data in Fig. 3 (a) are qualitatively 
identical with that in Fig. 2 (a): the correlated motion decays with the inter-particle distance R following the power 
law. As the area fraction n is increased, the curves shift to higher values and the decay rate of < D rr /r > and 

< Dgg /t > versus R becomes slow. 

Figure 3 (b) shows an anomalous long-rang correlation. For S3 particles, when n > 0.08 and R > 5d, the decay 
of < D rr /r > and < Dgg/r > becomes slow with increasing R and does not follow power law. At large R region, 

< D rr /r > and < Dgg/r > are weakly dependent on inter-particle distance R. For short inter-particle distance 
(R < 5d), the behaviors of < D rr /r > and < Dgg/r > are similar to that of SI and S2 particles. 



4 




FIG. 1: Correlation function < D rr /r > (solid symbols) and < Dgg/r > (open symbols) (scaled by Do) as a function of 
inter-particle distance R (scaled by d) with the area fraction n ranging from 0.01 to 0.05. The average result of these curves 
were drown in the inset, while the corresponding area fraction n was estimated as 0.03. The solid lines in the inset are the 
fitted line < D rr /r >= 0.48/.R a86±0 ' 02 and < Dge/r >= 0.28/ R 1A5±om . The samples are SI silica d = 1.57 ± 0.06/im. 




FIG. 2: (a) Correlation function < D rr /r > (solid symbols) and < Dgg/r > (open symbols) (scaled by Do) as a function of 
inter-particle distance R (scaled by d) with various values of area fraction n. Different symbols represent the data for different 
area fraction n. The data in the figure were fitted by < D rr /r >oc l/R Xr and < Dgg/r >oc 1/R X ° . (b) The fitted parameter 
A r and Xg as a function of the area fraction n for SI silica d — 1.57 ± 0.06/im. The samples are SI silica d — 1.57 ± 0.06/^m. 



III. DATA ANALYSIS 



Figure 2 (a), 3 (a) and (b) show that < D rr /r > and < Dgg/r > for various values of particle area fraction n are 
different. The exponent A r and Xg decrease monotonically with the area fraction n. Our results shared the similar 
trend of the work in Prasad et al. (Iol |. where the exponent is monotonically decreasing by adding HSA protein onto 
surface. In Prasad's system, the weight of surface viscosity rj s was modified by HSA. By scaling with ?y s , the whole 
dispersed curves could collapse onto a single master curve. In our data, D Q should be replaced by D s = ksT '/iirrjes, 
where r/ es is effective surface viscosity. The horizontal axis of R should be scaled by the ratio between the bulk viscosity 
and the surface viscocity. Following the scaling method of Prasad et al. [l(|, we replot the curves for different n by 
scaling the two sets of effective surface viscosity rj es and rj es for separation and correlated motion individually. The 
scaled correlation function and separation are follows: 



D" rr ,ee =< D rr fig/r > /2D S , S = Rf] w /-q e 



(2) 
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FIG. 3: Correlation function < D rr /r > (solid symbols) and < Dge/r > (open symbols) (scaled by Do) as a function of 
inter-particle distance R (scaled by d) with various values of area fraction n. Different symbols represent the data for different 
area fraction n. The samples are the following: (a) S2 silica d = 0.97 ± 0.05/im; (b) S3 silica d — 0.73 ± 0.04pim. 





FIG. 4: Master curve of the correlation function < D rr /r > (solid symbols) and < Dee/r > (open symbols) (scaled by 2D a ) 
as a function of reduced particle separation S = Rrj w /rj s with various values of area fraction n. Different symbols represent the 
data for different area fraction n. The data and symbols in Fig. 4 (a), (b) and (c) are same as that in Fig. 2 (a), 3 (a) and 
(b), respectively. The samples are the following: (a) SI silica d = 1.57 ± 0.06/im; (b) S2 silica d — 0.97 ± 0.05pim; (c) S3 silica 
d = 0.73 ± 0.04^m. 
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FIG. 5: The effective surface viscosity r) es versus rj es for particles SI, S2 and S3. The samples are the following: (a) SI silica 
d = 1.57 ± 0.06^m; (b) S2 silica d = 0.97 ± O.Odfim; (c) S3 silica d = 0.73 ± 0.04^m. The solid lines are the fitted curve 
v'L = K(r)' es ) 2/:i , where K = 0.88 ± 0.02, 0.84 ± 0.04 and 0.77 ± 0.03 for (a), (b) and (c), respectively. 



S = Rv w /Ves, (3) 

where D s = hsT/Anr)'^ and i] w is the viscosity of the water. We allow them to vary independently to obtain two 
independent measures of the effective surface viscosity r] es and r\ es as Ref.[10]. 

Figure 4 (a), (b) and (c) show the relationship of the scaled variables D" rr , D" gg to the scaled separation S with 
various values of n for samples SI, S2 and S3, respectively. The data and symbols of Fig. 4 (a), (b) and (c) are same 
as that in Fig. 2 (a), 3 (a) and (b). The curves with various values of n fall on a single master curve. The curves 
with higher values of n fall on the left side of the master curve. Figure 5 (a), (b) and (c) describe the scale factors 
ij es and ij es used to create the master curve for samples SI, S2 and S3 respectively. As shown in Fig. 5, rj es is clearly 
deviated from rj es for the same area fraction. And the discrepancy between r) es and rj es increases as the area fraction 
n is creased. The relationship of -q es to i] es for samples S2 and S3 is qualitatively same as that of particles SI. 
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FIG. 6: Master curve of the correlation function < D rr /r > (solid symbols) and < Deg/r > (open symbols) (scaled by 2D S ) 
as a function of reduced particle separation L with various values of area fraction n. Different symbols represent the data 
for different area fraction n. The data and symbols in Fig. 6 are same as that in Fig. 2 (a). The samples are SI silica 
d = 1.57 ± OMfim. 




FIG. 7: (a) The mergence of master curve for particles SI and S2. The scaled variables D rr =< D rr /r > /2D S (solid symbols) 
and Dge =< Dgo/r > /2D S (open symbols) as a function of reduced particle separation L for SI and S2 particles can collapse 
onto a single curve, (b) The mergence of master curve for particles SI and S3. The scaled variables D rr =< Drr/r > /2D S 
(solid symbols) and Dgg =< Dgg/r > /2D a (open symbols) as a function of reduced particle separation L for SI and S3 
particles can collapse onto a single curve. Different symbols represent the master curve of different particles. 



The results shown in Fig. 5 are hard to understand. The r] es and r] es should be identical, since they describe the 
same physics feature. However, experimental results show that the discrepancy between r] es and r\ es is extraordinarily 
large, e.g., it exceeds 50 percent for samples S2 while n = 0.2. In Rcf.[10], the deviation of rj es and r\ es was only within 
15 percent and their averaging values agree well with one-particle viscosity. In our case, the discrepancy is too large 
to fall in our experimental error range. Especially considering the dependence form of r\ es ~ (rj es ) 2 ^ 3 in caption of 
Fig. (5) , the exponent of 2/3 means that the dimension of ri es on rj es are not matching even, which is unreasonable in 
physics. Checking the above scaling process and comparing with r\ es ~ (r? es ) 2 ^ 3 , we think that the dimension scaling 
method rj w /ri es may be not necessary way. The fitted index 2/3 suggests that the following scale method maybe a 
reasonable choice: 

L=^(^) 3/ \ (4) 



and the correlation function D rr ^g =< D rr> gg/T > /2D S where D 8 — IcbT '/ '4irr] es . By this method, we could scale 
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n es ,i P < cP ^ m > 

FIG. 8: Two-particle surface viscosity r/ es ,2p versus one-particle surface viscosity r] e s,\ P for particles SI, S2 and S3. The samples 
are the following: (a) SI silica d = 1.57 ± 0.06^m; (b) S2 silica d = 0.97 ± 0.05^m; (c) S3 silica d = 0.73 ± 0.04^m. The slope 
of the solid lines is 1, indicating an equality between ?7 eSl 2p and r/ eay ip. 

the correlation function and separation onto a single master curve with one set of values of parameter rj es rather than 
two sets of parameter value. The obtained value of rj es by this method is called two-particle effective surface viscosity 

Figure 2(a) for samples SI was replotted in figure 6 with the new scaling method described in formula (4). The 
symbols are same as previous. As shown in Fig. 6, all the data sets fall on a single master curve. For samples S2, the 
behavior of D rr , Dgg is identical with that of samples SI and the data sets of the two kinds of particles can fall on 
the same master curve, as shown in Fig. 7 (a). For particles S3, in the case of short inter-particle distance R < 5d, 
the scaled correlated motion D rr , Dgg for different n can fall on a single master curve, whose behavior is similar to 
that of particles SI and S2. However, when R > 5d the deviation from the master curve appears, which is different 
from that of particles SI and S2. Figure 7 (b) shows that the data of scaled correlated motion D rr , Dgg of particle 
SI and S3 for the case of R < 5d can also collapse onto the same master curve. The identicalness of the master curves 
for different samples suggests that particle size dependence has also been eliminated. The superposition of the data 
sets for various values of n and different particles indicates that there is a general uniform relationship between cross 
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FIG. 9: Two-particle surface viscosity 7/ es ,2p as a function of area fraction n for particles SI, S2 and S3. Different symbols 
represent the data of different particles. The solid lines are the fitting curve rj es ,2p = %\J~^ ex P{—-^\J~-&i) + V^ro- The fitted 
values of the parameter r and k are presented in table I. 



correlation coefficient D rri Dgg and effective viscosity rj es . 

In the experiment, great care is taken to clean the interface to eliminate polymer, surfactant and contamination's 
effects. We have used Langmuir trough method to measure the changing of surface tension while the area of aqueous 
surface reduced. No detective impurity was found in the system. We believe that there is no heterogeneity in our 
experiment system. The two-particle effective surface viscosity r? e s,2p should be nearly equal to the one-particle 
effective surface viscosity r] eSi i p which is directly related to particles's self-diffusion. To compare r/ es ^ P and rj es> i p , we 
calculate one-particle effective surface viscosity r] eSy i p as follows. At Ref.[14] Chen et al., presented the short-time 
self-diffusion coefficient D s s of silica particles at an air- water interface via the formula < Ar 2 (r) >= 4D s s t for different 
particle concentration. The measured D s s /Dq can be fitted to a second-order polynomial, 



D S JD = a(l -0n- yn.% (5) 

where a, j3, 7 are the fitted parameters. At the n — > limit, the short-time self-diffusion coefficient D^(n) is directly 
related to the drag coefficient £ via the equation D s s {n = 0) = &bT/£. For particles at a distance z from the interface, 



t=(rjia)f(z/a,B), (6) 

where z is the distance between the sphere's north pole and the interface, i.e., z = when the sphere is in contact with 
the interface from below, and a is the radius of the sphere [l5|. The Boussinesq number B is defined as B =Jh/ {t] w o), 
where r) s is the shear viscosity of the interface and t] w is the viscosity of bulk of water. Fischer et al. [l5l. Ila| showed 
that f(z/a,B) = + k^B + o(B 2 ), which reads the effects of interface viscosity on the single-particle. According 
to the relations of short-time self-diffusion coefficient to the drag coefficient, the one-particle effective surface viscosity 
can be written as: 

Ves ' lp = (a(l-/3n- 7 n 2 ) ~ ^ ^ ' (?) 
The zeroth-order and first-order drag coefficient k^ and fc^ are given by 15]: 

k^^f^^), (8) 



fc (i) w I -41n(farctan(|))( I ^l 75 ) for z>0 
\ -4 In (| arctan(^t^)) for z<0. 
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The zeroth-order drag coefficient fcW w 15.97, 13.0 and 15.97 for samples SI, S2 and S3 respectively flij . From fe(°), 
we derived the values of fcW described by formula (9). Substituting k^°\ fcW and the measured values of a /3 and 7 
at Ref.[14] into formula (7), we derived the one-particle effective surface viscosity f] eSy i p for different area fraction n. 

Figure 8 (a), (b) and (c) exhibit the relationship of the one-particle effective surface viscosity r) es ,ip to the two- 
particle one r?es,2p used to create the master curve for particles SI, S2 and S3, respectively. As shown in Fig. 8, the one- 
and two-particle measurements agree well with each other. The results confirm that the non- heterogeneity argument 
about our experiment system and the scaling method described in formula (4) are reasonable. But at present we do 
not derive the scaling relationship of formula (4) theoretically, which deserves further investigation. We vaguely guess 
the index 3/2 may be related to the dimension feature of the effective surface viscosity, which dependents on both of 
the bulk and interface viscosities. 

Most interesting phenomenon in the experiment is that the exponent A r and Xg decrease when the colloidal area 
fraction n increases. Experiment results in Fig. 1 indicate that the cross correlation between particles in our experi- 
ment is stronger than that exhibited in the experiments of Prasad et al. [loj even in the case of low r/ s . By all above 
analysis, the decrease of X r and Xg value is correspond to the increase of the two-particle effective surface viscosity 
Ves,2p- Figure 9 shows the two-particle effective surface viscosity r) es ^p as a monotonous increasing function of the 
area fraction n for samples SI, S2 and S3. The values of n es- 2 P for particles S2 are extraordinary smaller than that of 
SI and S3. 



IV. DISCUSSION 



Silica spheres representing an important class of charged particles commonly used in colloidal science have anionic 
SiO~ groups on their surface [l7], HI]- These charged spheres at the air- water interface show a DLVO repulsive 
interaction [14j | . The interactions between silica spheres at air- water interface with low fraction n include hydrodynamic 
and Coulomb interactions. With the increase of n, the mean particle separation £ — d\Jn/ (An) decreases and the 
mean Coulomb interactions between particles increase. The Coulomb interaction between colloidal particles should 
serve as an effective surface viscoelastical role in our system and enhance the cross correlation of particles diffusion 
(shown as Fig. 6 and Fig. 9). Regarding the effective surface viscosity n eS:2p as a measurable scale of the energy-cost 
due to the deformation of interfacial particles configuration, the dependence of n eSt 2 P on n should share the feature of 
interaction potential among the silica sphere. 

Particles' interaction potential in our system can be well described by the screened Coulomb potential [l4| 



exp[-R/X D ] 
U C (R) oc - , (10) 

where R is the separation between the particles, and Xd is the screening length. Referring to this potential form and 
relationship between R and n, the effective surface viscosity n es ^ P could be written as: 



k [An d nr. 

Ves,2 P = -t\/ — exp{ \H-)+Vwr , (11) 

a V 7r ro V An 

where d is the particles' diameter and n is the area fraction. Two fitting parameters are k and ro- The parameter 
k is related to the strength of the influence of Coulomb interaction on the effective surface viscosity. The parameter 
ro is a intrinsic length of the interfacial system. Considering the definition of surface viscosity, ro maybe related to 
the length of the line that spheres contact with the interface. As an amplitude parameter, k should be related to 
the surface charge, immerse depth and diameter of the particle, though no exact expression of n couldn't be given at 
present. Here, we add an extra term of rj w ro as a correction to the theory of equation (7), since these silica spheres, 
especially samples SI and S3, experience a larger drag force than the value predicted by T.M.Fisher's theory flij . 
This term of rj w ro represents the surface viscosity that one particle experiences at the air-water interface when n — > 0. 

The solid lines in Fig. 9 are the fitted curves described by formula (11). The fitted values of parameters k and ro 
are given in table 1. The value of ro for samples SI is nearly 2 times larger than that of S3, and k close to 4 times in 
the error range. These ratio between numbers are consistent well with the size ratio between the particles, considering 
the underlying physical meaning of k and ro, viz. surface charge and immerse contact line. But for samples S2, the 
values in table give a different story: a much larger n vs. a much less ro. These two numbers are self consistent at 
least. More surface charge (^ k) will give a shorter contact line (^ ro), since high charged particles always immerse 
deeper. As mentioned above, the S2 particles were purchased from Bangs Laboratories, while SI and S3 particles 
from Duke Scientific. Both were synthesized in different ways fl4j . The difference in surface chemistry may introduce 
a different immersion ratio. 
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TABLE I: Particle sample used in the experiment and the fitted values of ro and k from the measured effective surface viscosity 



Sample/Manu. 


d(fim) 


ro(/xm) 




silica 1 (Sl)/Duke Sci. 


1.57±0.06 


1.13± 0.1 


510±84 


silica 2 (S2)/Bangs Lab 


0.97±0.05 


0.30± 0.05 


1500±800 


silica 3 (S3) /Duke Sci. 


0.73±0.04 


0.67± 0.1 


160±34 



At a glance, this S2 immerse-deeper picture is sort of conflict with the fact that the D s s /Dq value of S2 1.5) is 
higher than that of SI and S3 1.2) when n — > [l4j, viz. the particles S2 experience a relative less drag than SI 
and S3 on the air-water interface in the dilute limit. While n — > 0, the effective viscosity r\ e is composed of the bulk 
viscosity r\ w and interfacial viscosity r] es = rj w r . The short self-diffusion constant of particle could be written as: 



Dl(n -»• = = -, (12) 

Vw+Ves r] w ak^> + r) w r 

where a is the particle radius, and r is the character length. The deeper the particles immerse, the larger the bulk 
viscosity ij w contributes and the less the r] es . Base on our data, the increase of t] w might be overwhelmed by the 
decrease of rj es when S2 immerse deeper. 

The r in equation (12) should be same with the ro appeared in equation (11), since both describe the interfacial 
viscosity r] es when n — >■ 0. Regarding r as the length of interfacial particle contact line, we could calculate the value 
of r by substituting the fitting value of ro (in table I) and the value of D s s (in paper [lj]) into the formula (8) and 
(12). The calculation results show that r = 0.97 ± O.l^m, —0.87 ± 0.1/xm and 0.47 ± 0.1/j,m for samples SI, S2 and 
S3, respectively. The deviation between ro and r is less than 15 and 30 percent for sample SI and S3, respectively. 
However, the sign of the r is negative for S2. With a long shot guess, it maybe because that the particle S2 has a 
partially slip boundary condition due to its unique surface chemistry. Then the fitting value of ro is not the real length 
of particles contact line, but a viscosity-effect-equivalent one, which could much less than the real contact length. By 
this underestimated ro, we could have a overestimated rj w due to a deeper-immersion estimation of z/a in equation 
(8). In the end, we get a minus r for sample S2. But at present, we are unable to experimentally measure the exact 
length of the interfacial particle contact line independently. So we cannot say anything much more here. 

Traditionally, particle concentration dependence of the bulk viscosity (data points in Fig. 9) could be expressed by 
the polynomial OHO]. By analogy with the bulk case, we rewrite the effective surface viscosity rj es ^ P as follows: 



T]es,2p = TJesQ (1 ~ U) M , (13) 

where 77 es o and M are f ne fitting parameter. The parameter 7/ es o represents the effective surface viscosity that a sphere 
experiences. The area fraction dependence of r] eSi 2p is represented by /x. The data points in Fig. 9 could also fitted 
very well by equation (13). However, the fitted values of /i are very high and /i ~ 12 — > 15, indicating rj e s,2p has a 
super strong area fraction dependence. Such incredible high power exponent suggests that rj es ^p should exponentially 
increases with colloidal area fraction. We believe that the strong area fraction dependence of r\ es ,2p does arise from 
the Coulomb interaction between particles. 

Generally, figure 9 exhibits dependence of shear viscosity rj eS 2 P on particle concentration, which is directly related 
to cross-correlation of interfacial diffusion. With the increase of area fraction, the mean particle separation decreases 
and the electrostatic interaction between particles enhances. The enhancement of the electrostatic interaction causes 
the increase of correlated strength between particles. The increase of correlated strength corresponds to the increase 
of the effective surface viscosity r? es ,2p- We describe the electrostatic interaction between particles with the effective 
surface viscosity T\ ea> i v . There are hydrodynamic and electrostatic interaction between particles in our experimental 
system. We call this interaction effective hydrodynamic interaction. With the increase of n, the effective hydrodynamic 
interaction enhances mainly arising from the increase of Coulomb interaction between particles. H. Diamant et al. 
0)HI| had demonstrated that many-body effects do not influence the hydrodynamic coupling at large distance in quasi- 
two-dimension (Q2D) colloidal system for momentum diffusion does not contribute to the large-distance coupling. 
We also observed that the hydrodynamic interaction between hard sphere is independent of particle concentration in 
the two-dimension (2D) system where the momentum leaks to the third dimension (22j. 

Another information given by Fig. 9 is size dependence of rj es ^p- First, we study the behavior of the shear viscosity 
Ves,2p for particles SI and S3 both which were purchased from Duke Scientific with identical surface chemistry. It is 
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shown that for the same area fraction the shear viscosity rj eSt 2 P of small particles is slightly larger than that of large 
particles, i.e., small particles experience a larger drag than large particles. One possibility is that smaller particles 
are separated by a smaller distance (I = dy^ir/ (4n)) as compared to the larger particles at the same area fraction n, 
and they may feel more effective hydrodynamic interactions at air-water interface than the larger particles do (23l - [25| . 
Compared with particles SI and S3 ,the S2 particles has very small shear viscosity rj es ^p, as mentioned above, which 
maybe come from the different surface chemistry. 

Figure 3 (b) shows an anomalous long-rang correlation. While R > 5d, the decay of < D rr /r > and < Dgg/r > 
bates with increasing R and does not follow power law. At large inter-particle region, correlated motion are weakly 
dependent on R and are roughly logarithmic (too narrow data range to tell, indeed). The deviation of correlated 
motion from the master curve in the case of R > 5d indicates another rules governed the system or something else? 
At present, we don't have an answer yet. During the experiment, we use 10-100 image sequences, each contains 
100 images, to calculate D rr and Dgg, and the result is further averaged over repeated runs (10—20 runs). This 
corresponds to an average over 10 6 - 10 7 particles, ensuring that the statistical averaging is adequate. We conduct the 
measurements at different times with separately prepared air-water interface. The identical experiment equipments 
and measurement methods are used for all the three particle samples. Thus, this anomalous long-rang interaction 
sounds unlikely caused by inaccurate measurements. But this anomalous correlation occurs only at smallest colloidal 
S3 system. The diffusion constants of S3 particles is relative larger than two others samples due to the smaller 
size. An improved microscopy system with even higher time and spacial resolution might be necessary to clarify this 
phenomenon. 

V. SUMMARY 

In this paper, we systematically investigated the spatial correlated motion of weakly charged silica spheres at a 
pure air-water interface. Three kinds of silica particles are used in the experiment: silica 1 (SI) with the diameter 
d = 1.57(im purchased from Duke Scientific, silica 2 (S2) with the diameter d = 0.97/im purchased from Bangs 
Laboratories and silica 3 (S3) with the diameter d — 0.73/im purchased from Duke Scientific. Optical microscopy 
and multi-particle tracking are used to measure D rr /r and Dgg/r as a function of the inter-particle distance R with 
various values of particle concentration for particles SI, S2 and S3. Correlation function D rr /r and Dgg/r decay 
with inter-particle distance R as iy^°- 86 ± - 02 ancl l/R 1A5±0 - 03 : respectively, for particles SI with n = 0.03. With 
the increase of n, the decay of correlated motion with R becomes slow (Fig. 1 and 2). We attribute it as that the 
cross correlation is enhanced by the Coulomb interaction between particles. We describe the electrostatic interaction 
between particles at the interface with the effective surface viscosity r\ ea . 

Following the scaling method given by Prasad et al. 10], we scale the curves with various values of n onto a single 
master curve by T) es and r\ es . The measurements results show that r\ es is clearly deviated from ij es for a certain area 
fraction (Fig. 5). The discrepancy between r] es and rj es even exceeds 50 percent in the case of large n. The dependence 
of rj es on r] es suggests that we should rescale the separation according to formula (4) . With this new scaling method 
the curves with different area fraction can be scaled onto a single master curve by one set of parameter values of 
two-particle effective viscosity r\ ea ^. P (Fig. 6). The measurements of effective surface viscosity rj es ^ P and rj es ^i p agree 
well with each other (Fig. 8). The behaviors of cross correlation for particles SI, S2 and S3 are similar and their 
master curve can also fall on a single curve except the large R region of S3 (Fig. 7). The experimental results indicate 
that our new scaling method is reasonable and could be used in similar experiment system. But presently we do 
not deduce theoretically the scaling relationship described in formula (4), which deserves further investigation, and 
also cannot understand the anomalous long-rang correlation exhibited by particles S3 at large inter-particle distance 
region (Fig. 3 (b)). 

The effective surface viscosity rj es ^ p as a function of n for samples SI, S2 and S3 are presented (Fig. 9). For samples 
SI and S3 with identical surface chemistry, rj es ^p of small particles is slightly larger than that of large one for the 
same area fraction. However, the viscosity rj es ^ p of particles S2 is extraordinary small compared with that of SI and 
S3. The difference of surface chemistry characteristic and immersion depth between particles SI, S3 and S2 may be 
the two possibilities causing the discrepancy of the viscosity r\ es ^. P - The experimental data about the effective surface 
viscosity rj eSt 2 P presented in the paper are helpful for the related theoretical investigation. Referring to the form of 
the interaction potential between particles, we present a exponential fitting formula (equation (11)), which is more 
physically reasonable than the traditional polynomial form, for the viscosity r] es _2 P - 

The decay of the correlated motion with R for lowest n in this paper is slower than that of Ref.[10] with low surface 
viscosity. We argue that the high sensitivity of electrostatic interaction between particles on area fraction results in 
this deviation. With the increase of n, the electrostatic interaction enhances and the slope of correlated motion vs. 
R decreases, which confirm our argument. We increase the Coulomb interaction between particle by increasing the 
area fraction n. Though we have known that many-body effects would not influence the hydrodynamic interaction 
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of particles, the method varying Coulomb interaction between particles in the paper is indirect. We will change the 
surface charge of particles directly to investigate the effects of Coulomb interaction on correlated motion in the next 
experiment to further confirm our argument. 
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